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Based on differences of generalized Renyi entropies nontrivial constraints on the shape of the 
distribution function of broadly distributed observables are derived introducing a new parameter 
in order to quantify the deviation from lognormality. As a test example the properties of the 
two-measure random Cantor set are calculated exactly and finally using the results of numerical 
simulations the distribution of the eigenvector components calculated in the critical region of the 
lowest Landau-band is analyzed. 
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Numerical simulations of disordered critical systems 
always face the problem of dealing with system sizes 
smaller then the relevant length scale therefore highly 
complex fluctuations appear. Consequently expectation 
values of observables do not represent the full picture for 
the description of the dynamics and physics involved, in- 
stead the whole distribution function of these observables 
should be studied. This problem has become the central 
issue of theoretical research over the past decade . 

Multifractal analysis has now become a standard 
tool to handle the aforementioned fluctuations aris- 
ing near criticality. Especially the disorder induced 
localization-delocalization (LD) transition is an impor- 
tant and widely studied example of such critical sys- 
tems d). 

The general method of multifractal analysis starts with 
the calculation of the measure Q{i) of the observable over 
boxes of linear size in a system of length L divided 
into M pieces of these boxes. In the thermodynamic 
limit of i/L — > the qth moment of this distribution 
scales with a nontrivial exponent Dq that is, in general, 
strongly g-dependent. These moments are often called 
as generalized inverse participation numbers (IPN). The 
quantity (q — l)Dq, also known as the mass exponent r^, 
is a monotonous, nonlinear function of q. 

Another way of describing these fluctuations is using 
the f{a) function, the Legendre-transform of Tq. The 
/(a) is a smooth, single-humped, concave function. It is 
the distribution of the exponents aq — dr/dq that charac- 
terize the fiq{£) ^ {£/L)°''' scaling of the box-observables 
where fj,q {£) is the total measure of the qth power of the 
observables Q over the boxes with linear length £. The 
f{oiq) gives the scaling of the number of boxes with mea- 
sure (IjLy^i in the interval (a, a -I- dct) 

p(t,a)=5(t,a)e*/(") (1) 

where we have introduced the variable t = — ln(£/L). 
The above distribution describes the multifractal in the 



sense that lnp(i, a)//: — > /(a) if i ^ oo. The /(a) func- 
tion obeys certain general rules: /'(a^) — q, where prime 
denotes derivation with respect to a. For some specific 
values of q 

f'iao) = 0, /(ao) = Do (2a) 
/'(ai) = l, /(ai) = ai=L»i (2b) 

where Do is the Hausdorff-dimcnsion of the support of 
the observable Q. Furthermore the support of the f{a) 
spectrum is a finite interval [ 1, and at the two 

limiting values the derivative of the f{a) is infinite. The 
importance of the study of the function /(a) is enhanced 
due to its connection to the properties of the local distri- 
bution function of the Q variables as 

n(i,lnQ)-exp[i(i?o-/(a))], a = lnQ/t (3) 

The shape of the /(a) curve is close to a parabola 
which means that the distribution Tl{t, In Q) is a Gaus- 
sian and the n(i, Q) is lognormal. Therefore any devi- 
ation from the parabolic approximation (PA) shows up 
as a deviation of 11 (t, Q) from lognormality. Numerical 
experiments show that there is a substantial deviation 
of the distribution function (||) from the simple Gaus- 
sian especially for large values of InQ, hence it is even 
asymmetrical (see Fig. |^). 

The calculation of /(a) and Dq functions in numerical 
experiments suffers from several difficulties. First this 
analysis is suitable only if £ ^ L < ^ is satisfied, where ^ 
is the localization length and still the box-length should 
be larger that any microscopic length scale of the system. 
Second, higher moments of the distribution, \q\ > 3 are 
usually calculated with low precision and the PA is usu- 
ally satisfactory. The PA is fixed by the values of ao and 
Dq. In spite of the simplicity of the PA it breaks down 
at values of q where the monotonicity of the Tq is lost. 
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FIG. 1. The joint probability distribution function of the 
local amplitudes Q = [i/jp of the eigenstates taken in the 
critical region in a Quantum Hall System. The dashed line 
represents the parabolic approximation, the solid one is the 
improvement given in the text. 

Apart from the results of numerical simulations, re- 
cently it has been shown analytically that the moments 
of the distribution function of the eigenstate components 
in a two dimensional disordered system shows multifrac- 
tality however, only the PA has been derived. In 
another work the multifractal properties of a relativistic 
fermion moving in a two dimensional random vector po- 
tential field has been derived [0 without the restrictions 
of the PA. Therefore it seems to be of growing importance 
of deriving as much information from numerical experi- 
ments as possible. 

In this Letter we show that the deviations from the 
simple lognormal distribution can be measured with the 
application of specific Renyi-entropies [ pl| . These Renyi 
entropies are proportional to the variable t only for 
strictly self-similar observables in the thermodynamic 
limit, e.g. deterministic multifractals such as the two- 
measure Cantor set, however, they are connected to some 
specific values of the scaling dimensions of the general- 
ized IPNs and hence impose new conditions on the shape 
of the f{a) function. 

Hereby we use two differences of Renyi entropies [0 . 
The first one describes the deviation of the generalized 
IPN In P2 of the distribution Eq. (|l|) from the Shannon 
entropy H (a special Renyi entropy) and the other one is 
the normalized value of In P2 



The calculation of the quantities S2 and r7i2 involve in- 
tegration over the measure of the distribution in question. 
In our case, we have to perform integration of Eq. (|l|) over 
the allowed range of the a values [q;_oo, ct-i-00] 

Pi\e) - / f°'p{l,a)da (5a) 
H{£) Inl arp{e,a)da (5b) 



Furthermore, being interested only in the small £, i.e. 
large t limit, we can apply an asymptotic method due 
originally to Laplace [|6| and obtain an expansion of t of 
the form |17 



S^it) = (i^l -Z?2)t + ^S2,fci^'= 

00 

m2(i) = {D- D2)t + J2 m2M~^ 



(6a) 
(6b) 



fe=0 



From Eqs. (|6|) it is clear that their ratio tends to a con- 
stant as t — > 00 



S2{t) 
t->oo 1712(1) 



D1-D2 
Dn-D2 



(7) 



This is the key result of our Letter. The Di and D2 val- 
ues in Eqs. (||,0) are related to the /(a) function, since 
(|^) ai = /(ai) = Di and D2 = 2a2 - /(a2). On the 
other hand one can calculate the quantities 82/^12 di- 
rectly using the observable Q at any length scale. The 
latter value imposes a further condition on the f{a) 
through the dimensions Di and D2 using relation Eq. (|^): 
/(q!2) = 2a2 ~ {cti ~ J2D0) / {1 — 72). Note that for a log- 
normal distribution 72 = 1/2 independently of the details 
of the distribution as it is also 72 = 1/2 for a parabolic 
/(a), hence a deviation from the value 1/2 is a sign of 
non-lognormal distribution and non-parabolic f{a) [^ . 
Furthermore, the denominator of Eq. (|^) is the quantity 
that has been recently connected to the anomalous dif- 
fusion exponent by Refs. j|] and [p^ . 

As an illustration first let us consider the two-measure 
Cantor set representing a binary process . In this con- 
struction one divides the [0, 1] interval into two equal 
parts receiving different, p and 1 — p measures. While it- 
erating this procedure for each subhalf of the interval we 
get a multifractal in Dq = 1. This is a strictly selfsimilar 
multifractal if we keep the < p < 1/2 value fixed, be- 



S2 = H - In P2 and 777,2 = In A/ - lnP2 (4) cause for the nth stage = "-'5'2^' and 



,(1) 



Generalizations for higher moments are also possi- 
ble p^ , however, numerical simulations for such cases 
are presently not too reliable. The above two parame- 
ters have already been successfully applied in a number 
of systems [|l^,Q^ for the shape-analysis of the com- 
plex distribution function of eigenvector components and 
energy spacing distributions. 



where S'2^^ and 



,(1) 



are the values for the initial config- 
uration. Hence all the S2,n and the 7772, n in Eq- @ are 
zero, and the 72 quantity is independent of n and is a 
function of p 1 17|. 

Next we consider a generalization of this Cantor set 
by choosing at every stage a different value for p„ then 
we obtain a non-selfsimiliar multifractal the random bi- 
nary Cantor set (RBCS). It is straightforward to derive 
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the f{a) function once the distribution function of the 
Pn sequence is further specified. Taking a uniform distri- 
bution over the interval [1/2 — 6,1/2 + 6] we have found 
that 72 (<5) is a smooth function of < 6 < 1/2. Using 
this simple mathematical construction we will show that 
the PA can be easily improved with the application of 
parameter 72 Eq. (^. As an example we took 6 = 1/3. 
The PA in Fig. |^ is denoted as a dotted line while the 
exact relation is represented by the solid symbols. We 
can see that adding a fourth order term (a — ao )^ to the 
PA fulfilling both Eqs. (|) and (0) in /(a) gives the solid 
line that is clearly a much better approximation than the 
simple PA. 




FIG. 2. The distribution of the random binary Cantor 
set (RBCS) with S = 1/3. The dashed line stands for the 
parabolic approximation and the solid one is the fourth order 
one using parameter 72 (see. Eq. (0)). The symbols represent 
the exact relation. 

Finally we will apply our analysis to the widely studied 
LD transition that takes place in a two dimensional disor- 
dered system subject to a perpendicular strong magnetic 
field: the quantum Hall system [|6|jl9[|. This LD tran- 
sition is responsible to the integer quantum Hall effect. 
It is described by the critical index v of the localization 
length and the f{a) spectrum, that is believed to be uni- 
versal 1^, however, up to now the focus has been set on 
the PA only and the value of ao fixing the position of 
the maximum of it. The value of ao describes the scaling 
behavior of the typical local electron density. We show 
how the deviation from the lognormal distribution can 
be measured using Eq. (0) and propose an optimal mod- 
ification of both the f{a) and the distribution function 
H(i, InQ). For that purpose we have performed calcu- 
lations of two dimensional systems of linear size of 200 
magnetic lengths. We have obtained 134 eigenfunctions 
and calculated the overall joint probability distribution 
function of the local amplitudes Q = In Fig. |l| we 

have plotted the logarithm of the histogram of InQ as 
a function of InQ/t. On the same figure the PA is also 
presented. The discrepancy is clear especially for larger 



values of Q. On the same figure we have given a possible 
improvement of the analytical /(a) that shows a much 
better resemblance to the numerical data. 
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FIG. 3. Scaling dependence of parameters m2(t)/t and 
S2{t)/t as a function of 1/t. Solid symbols stand for the val- 
ues obtained in the Quantum Hall System and continuous 
curves show several [K/L] order Fade approximants: dashed 
line [0/2], dotted line [0/3] and solid line [0/4]. 

The improvement in Fig. |l| was achieved similarly as 
for the case of the RBCS. For that purpose we have cal- 
culated the m2{t)/t and S2{t)/t values and plotted as a 
function of 1/t. Unfortunatelly the low system size is 
a severe limitation on our problem, since it is desired 
to obtain the values of m2{t)/t and S2{t)/t for t oo. 
However, it is not possible to go beyond t^ax = In L 
and current computation possibilities impose a limit of 
tmax ~ 6.0. In order to obtain an acceptable value for 
the values TO2 (t) /t and 5*2 (t) /i as 1/i — > oo we have per- 
formed an [K/L] order Fade approximation j2^ for both 
quantities. Besides finite size scaling this type of approx- 
imation is ready to give an estimate of the coefficients in 
the expansions Eq. ^ . Not having data obtained for dif- 
ferent system sizes at hand we can show that the Fade ap- 
proximation indeed gives acceptable results, as well. We 
assumed that the [K/L] order approximation should have 
no singularity in the given interval and has to fall to zero 
for 1/t ^ oo since on the scale £ ^ L the distribution is 
homogeneous, thus reducing to the possibilities to K < L 
and K = 0. In Fig. ^ we have plotted the calculated 
quantities together with several Fade approximants. The 
limit 1/i — > gives the approximate value for the param- 
eter 72. This way we obtained 0.445 for [0/2] 0.439 for 
[0/3] and 0.436 for [0/4]. Therefore we set the value of 
72 ~ 0.44. On the other hand, using standard multi- 
fractal analysis, our data gives ao =< InQ > /t ~ 2.24, 
Di = 1.78 D2 = 1.61 that are all well in the generally 
accepted range [p|p9|]. The latter values yield 72 ~ 0.44 
again. 

Having ao and 72 at hand we may extend the func- 
tional form of f{a) beyond the PA in order to get a bet- 
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ter aggreement with the numerical distribution function 
(see Eq. (^)). For the improvement of the functional form 
of the f{a) entering in Eq. (^) there should be a num- 
ber of possibilities. We used the addition of the term 
{a — ao)^ and (a — ao)^ to the PA fixing their coefficients 
requiring the constraints Eqs. (|^) and within an er- 
ror of 10~^. Therefore the n(x) function seems to be 
modified as U{x) — Upa{x) exp[— <(c3a;'^ -I- c^x^)] where 
X — InQ/t — ckQ, C3 « —0.03 and C5 « 0.12. 

In this Letter we have introduced the differences of 
Renyi entropies in order to describe the behavior of any 
multifractal as a function of t = — ln(£/L). These pa- 
rameters enable to look beyond the parabolic approxi- 
mation and obtain further constraints on the shape of 
the /(a) function as well as the distribution function of 
the multifractal quantity. In order to get insight in the t 
dependence we have applied Laplace's method and tried 
the results on exactly known deterministic and random 
multifractals. In the former case the self similarity prop- 
erty is also obtained while the latter case was a good test 
how the improvement from PA can be done. 

Finally we have analyzed the joint distribution func- 
tion of the local amplitudes of wave functions calculated 
in the critical regime of the Quantum Hall System. We 
have applied the Pade approximation in order to mimic 
the thermodynamic limit. The results show a possi- 
ble tail of the distribution function that is irregular for 
smaller amplitudes. 
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